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Abstract 

The aim of this paper is first to review the derivation of a model de- 
scribing the propagation of an optical wave in a photorefractive medium 
and to present various mathematical results on this model: Cauchy prob- 
lem, solitary waves. 

Keywords: Photorefractive media; Cauchy problem; solitary waves. 
AMS Subject Classification: 35A05, 35A07, 35Q55, 35Q60, 78A60 
I 

1 Introduction 

A modification of the refraction index in LiNb0 3 or LiTa0 3 crystals has been 
observed in the 1960s and first considered as a drawback. This photo-induced 
variation of the index is called the photorefractive effect and occurs in any 
electro-optical or photoconductive crystal. Applications have been found in the 
1970s-1980s to real-time signal processing, phase conjugation, or amplification 
of beams or images. 

In this paper we are interested in deriving a not-too-simple but tractable 
mathematical model for the propagation of light in such materials. Solitonic 
propagation is one of our concern but we focus here on initial value problems. 
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A very complete review of solitonic propagation in photorefractive media may 
be found in [6] . Our derivation follows the same guidelines as theirs but point 
out the different approximations made for future mathematical studies. 

The outline of the paper is the following. In Sec. [2] we first derive the 
Kukhtarev model for the material and then couple it to a wave propagation 
model for light to obtain a complete set of equations. A ID model is obtained 
keeping only one of the two transverse space variables. This is a saturated 
nonlinear Schrodinger equation, the mathematical theory of which is addressed 
in Sec. [3]in arbitrary dimension. Section|4]is devoted to the study of the full 2D 
model with emphasis on the Cauchy problem and the solitary wave solutions. 

2 Derivation of the Model 

2.1 The photorefractive effect 

The propagation of an optical wave in insulating or semi-insulating electro- 
optical crystals induces a charge transfer. The new distribution of charges in- 
duces in turn an electric field which produces a variation of the refraction index. 
The main characteristics of this effect are the following: (1) Sensibility to energy 
(and not to the electric field), (2) Nonlocal effect (charge distributions and the 
electric field are not located at the same position), (3) Inertia (charges need 
a certain time to move), (4) Memory and reversibility (in the dark the space 
charge, and therefore the index variation, is persistent but an uniform light 
redistributes uniformly all charges — this yields applications to holography). 

The sensibility to energy reminds us of Kerr media yielding the classical cu- 
bic nonlinear Schrodinger (NLS) equation. The nonlocal effects will of course 
complicate the mathematical analysis compared to NLS equations, but the gen- 
eral ideas will be the same. In our final model, inertial effects will be neglected 
since time is removed from the material equations. Memory and reversibility 
effects involve ion displacement in materials like Bi 2 Te0 5 , which we will not 
take into account in the present study. 

2.2 The Kukhtarev model 

The physical modeling of the photorefractive effect assumes that charges are 
trapped in impurities or defaults of the crystal mesh. We chose here to derive 
the model only in the case when charges are electrons. Some materials like 
semi-conductors necessitate to model both electrons and holes. Therefore we 
restrict our study to insulating media. 

2.2.1 Charge equation 

Electrons come from donor sites with density Nd. This density is supposed to 
be much greater than that of the acceptor sites (impurities) which we denote by 
Na- The density of donor sites which are indeed ionized is iVjJ and we of course 
have < Na <C Nd- Local neutrality, i.e. no electrons in the conduction 
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band, corresponds to the relation — Na- The total charge is given by 

p = e (N+ -N A - n c ), (2.1) 
where e is the electron charge and n e the electron density. 

2.2.2 Evolution of ionized donor sites 

Photoionization and recombination affect the density of ionized donor sites. 
Photoionization is proportional to the density of not ionized donor sites (Ad — 
N^). In the dark it is proportional to a thermal excitation rate (3 but is also 
sensitive to light intensity I em with a photoexcitation coefficient s. Recombi- 
nation is proportional to the density of electrons and occurs over a time scale 
r = l/(j T N^) which does not depend on n e if the excitation rate is low, there- 
fore the total evolution of ionized donor sites is 

d t N+ = 09+ sI cm )(N D - N+) - 7r n e A+. (2.2) 

2.2.3 Charge transport 

Now the main point is to describe the three phenomena which contribute to 
the charge transport or current density. The first phenomenon is isotropic and 
is due to thermal diffusion. It is proportional to the gradient of the electron 
density. The electron mobility is denoted by /i, T is the temperature and k# 
the Boltzmann constant. The second phenomenon is drift and is collinear to 
the electric field E tot . Finallly, the photovoltaic effect is collinear to the optical 
axis c and proportional to the non-ionized donor density and the field intensity 
with a photovoltaic coefficient /3 p h- The total current density is therefore 

J = e/m e E tot + [ik B TS7n c + f3 ph (N D - N+)cI cm . (2.3) 

2.2.4 Closure of the model 

The closure of the model is first based on charge conservation and the Poisson 
equation: 

ftp + V-J = 0, (2.4) 

V • (e eE sc ) = p. (2.5) 

The crystal is anisotropic and this is accounted for in the relative permittivity e 
which is a tensor. A careful analysis of the different electric fields has to be 
done. In Poisson equation (|2.5[) . E sc is the space charge field which is induced 
by the charge density. The total field E to t only occurs in the equations through 
its gradient (Eqs. (12. 4ft and (|2.3|) ). Two fields are constant and disappear in 
the final equations: the photovoltaic field E p h = /3 p h7 r AAc/e/is = -EphC, and 
an external field E ox t which is often applied in one transverse direction on the 
faces of the crystal. A last contribution to the total field is E, connected to the 
light which propagates in the crystal and its description is given in Sec. 12.31 

The set of five equations (|2. 1[) — (I2.5[) is called the Kukhtarev model and was 
first given in [9]. 
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2.3 Propagation of the light wave in the crystal 

We have already introduced the relative permittivity tensor £ which plays a role 
in the description of the propagation of a light wave in the crystal via the wave 
equation: 

2 {eE) - c 2 V 2 E = 0. 

In a non-centrosymmetric crystal the preponderant nonlinear effect is the Pock- 
els effect which yields the following E-dependence for the permittivity tensor: 

e(E) = e(0) - e(f^E)e = n 2 - e(f^E)e, 

where r is the linear electro-optic tensor and n the mean refraction index. We 
now suppose that E is a space perturbation of a plane wave (paraxial approxi- 
mation) of frequency to, wave vector k and polarization e: 

E(t, x) = A(x) exp(i(ujt - k ■ x))e. 

Such a wave with polarization e only "sees" a part of tensor e(E), or equivalently 
a variation 6n of the refraction index n: 

Sn = fee re ee*lE. 

2n L 1 

Now we can write an equation for the amplitude A which takes into account the 
dispersion relation c 2 |k| 2 = n 2 uj 2 and the slowly varying envelope approximation 
in the k direction. We denote by Vj_ the gradient in the perpendicular directions 
to k and 

V 2 ^k- V + 2|k| 2 — 

n 



A(x)e = 0. (2.6) 



Of course, we can consider the superposition of such waves to describe for ex- 
ample pump and probe experiments. 

The system is now closed but it is impossible to solve Eqs. (|2.ip - (|2.6[) . We 
have to simplify them taking into account characteristic scales. Our description 
follows (or more precisely makes explicit the assumptions in [T^]) and is purely 
formal. The rigorous justification is certainly difficult and should include the 
approximations made in Sec. 12.31 



2.4 Characteristic values 

We first want to define a characteristic electron intensity no by considering 
uniform solutions in space and time. Equations (|2.ip . (|2.2[) and (|2.5[) yield 

p = e(N+ — N\ — n c ) = and = (ft + sI){N D - N+) - lr n c N+ . 

With a characteristic intensity Iq, neglecting [3 and assuming n c -C A^a, we have 
n = sI (Nt) - N A )/j T N A . 

There are three characteristic times: (1) the characteristic lifetime of an 
electron (in the dark) r c = I/^-Na, (2) the characteristic evolution time of 
ionized donors = I/j^hq (and a consequence of no «C N A is r c <C Ta), (3) 
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the characteristic relaxation time of the electric field to — eo£ c / e.p,na, where e c 
is the characteristic value of e along the c direction. It is obtained combining 
Eqs. (|2.3[) and (|2.5p assuming there is only drift. If a timescale has to be kept, 
it is £o j but we do not detail this point since we neglect time-dependence in the 
final equations. 

The Debye length Ld is the characteristic value of the field space variation. 
It is determined together with the characteristic field Eq . The Poisson equation 
(|2.5p yields L B = £q£ c Eq/ eN A . If drift and isotropic diffusion have the same 
order, E = k^T/eL-a and therefore 

fkBTeo£c\ 1/2 _ f k B TN A \ 1/2 k B TN A 

L ° = — 2^7 — , Eo= \ > and I = . 

V e z N A ) \ e e c J s e c 

2.5 The Zozulya— Anderson model 

Zozulya-Anderson model 16J is obtained using the above characteristic values 
and for a specific material (LiNb0 3 ) which imposes certain symmetries. The 
adiabatic assumption allows to get rid of the time-dependence and an asymptotic 
formal analysis which accounts for the very large donors density Nn ends the 
derivation. 

Dimcnsionless equations are obtained using no and N A for electron and ion 
densities respectively, Iq for intensities, Eq for fields and e c for the permittivity 
tensor. Coefficient (3 is normalized as a dark intensity = (3/sIq. We keep all 
the other notations but they now denote the normalized variables. The total 
intensity is I = I em + id- We assume that the space charge field E sc derives 
from a potential: LdV</? = — E sc . In the adiabatic assumption matter equations 
reduce to 

1-N+N A /N 
1-N A /N D -" ciV °' 



LdV • < n c E to t + iDVn, 



1-N+N A /Ny \ 

+ E ph 7 em 1 _ Na/Nu )-0, 



-LlV ■ (eVip) = N+-1- ^-n D . 

In LiNb0 3 , N A /N D - 1(T 3 and n /N D - 1CT 6 and we neglect them. Finally, 
we make different assumptions on the fields: first the beam is not too thin, the 
photogalvanic and the external applied fields are not too large and therefore we 
may neglect — L^V • (eVip); second the propagation field amplitude is relatively 
small and we assimilate E to t and E sc . This implies n c = I and = 1 and we 
have only one matter equation, namely 

VI • X7tp + IX7 2 (p - V 2 / - k D E ph c ■ VI = 0, 

where /cd = V-^D- To obtain a "simpler" equation, in physics papers the 
variable U — ip — In / is often used. This variable seems however to lack physical 
meaning. 



5 



The final matter equation is 



VU ■ V<£ + V 2 C/ - k D E ph c • VI = 



(2.7) 



We now fix different space directions. Propagation is supposed to take place 
in the z-direction and k = ke z . The two transverse directions are therefore 
x and y. The direction is chosen as both c and e. If an external field is 
applied, it will be along e x as well. In the matter equation (|2.7J) . the quantity 



c • V/ simply reads d x I. In LiNbO„ 



is responsible for the change of 



refractive index (it is r xxy in some other materials) and we approximate e by 
n 2 in the expression for 6n which becomes Sn — ^n 3 rEQLj^d x (p. Together with 



Eq. (|2.6p the envelope equation now reads 



A(x) 



k 

-i—n 2 rEQL^d x ipA(x.) 



The last step is to have dimensionless space variables. We set a = ^n 2 rEo 
which has the dimension of the inverse of a space variable. We denote z' = \a\z, 
{x',y') = y/k\a\{x,y), A' = A/^h, <p' = vW*>/*d and U' = ^k\^\U/k B . 
The last approximations are now U' = tp' , k ^> \a\ and i? p h ~ Eq, and omitting 
primes: 



A = 



-iAd x ip, 



Vl^ + V ± ln(l + |A| 2 )-V ± ^ = d x ln(l + \A\ 2 ). 

These equations are usually referred to as a model derived in [16j but only 
seeds of these equations are derived there usually including many other terms 
and especially time derivatives. 

In the wide literature devoted to photorefractive media, many equations are 
written which resemble those above but with different choices of asymptotic 
approximations. In particular numerical results are very often obtained keeping 
the time in the matter equations (see [T3] or [TIj]). 



2.6 Mathematical setting 

If we look at a wider class of materials we may have different signs for the 
nonlinearity (in reference to the cubic nonlinear Schrodinger equation, the case 
a = 1 is classically called the focusing case, and a = — 1 the defocusing case). 
Besides mathematicians are more accustomed to use t as the evolution variable. 
We will therefore consider the system 

J id t A + AA = -aAd x (p, a = ±1, 

\ A^ + Vln(l + |A| 2 )-V^ = cUn(l + |A| 2 ), [Z ' S) 

where A = d 2 + d 2 or A = d 2 . 

These expressions with logarithms are widely used in the physics literature, 
maybe because they are the starting point of solitonic studies and logarithms 
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appear naturally in the expression of solitary waves (see Sec. 13. 2p . This form 
is however cumbersome to handle for the mathematical analysis, and it is much 
more convenient to cast (12.81) as 



( id t A + AA = -aAd x ip, 

{ div ((1 + L4| 2 rvV) = d x (\A\% 

which is closer to the original Kukhtarev equations. 

We have seen that the main effects take place in the t- (propagation) and 
the x-directions (drift, anisotropic diffusion, external field, polarization). It is 
therefore natural to study the equations with no dependence in the y variable. 
In the one-dimensional case, we infer immediately from the last equation in 
System (1231) that (1 + \A\ 2 )d x ip = \A\ 2 - C(t) where the constant C(t) is given 
by the boundary conditions. If no external field is applied C(t) = 0. This is 
the case for bright solitary waves (see [11]). In the case of dark solitary waves 
C(t) = \im x ^± OQ \A\ 2 (see [12]). which does not depend on t either. In both 
cases, System (|2.9p reduces to the saturated NLS equation 

l d t A + dlA=-a ^ A [-]^ ? A. (2.10) 

In the sequel we will mainly consider the case when A^ = and show that, 
in some sense, the dynamics of (|2.9[) is similar to that of (|2. 10[) which we will 
recall in Sec. [3] 

In the two-dimensional case (|2.9[) can be viewed as a saturated version of 
a Davey-Stewartson system. Namely, replacing 1 + \A\ 2 by 1 in the L.H.S. of 
(JHU) yields 

f idtA + AA = -aAd x tp, 
{ Aip = d x (\A\ 2 ), 

which is the Davey-Stewartson system of the elliptic-elliptic type (see [7]). 



3 The Saturated NLS Equation 

We review here some mathematical facts, more or less known, on the saturated 
NLS equation 

\A\ 2 A 

ldtA + AA = -a^ w , a = ±l, (3U) 
A(x,0) = A (x), 

where A = A(x, t) and x G R d . We have derived this equation for d = 1, but 
give here results for a general d. This equation is also derived in other contexts, 
for example the propagation of a laser beam in gas vapors [14] . 

3.1 The Cauchy problem 

The Cauchy problem (|3.1ip can be solved in L 2 and in the energy space H 1 . 
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Theorem 1 (i) Let Ao £ L 2 (R d ). Then there exists a unique solution A £ 
l;L 2 (R d )) of (|3.11[) which satisfies furthermore 



|A(x,t)| 2 dx 



|A,(x)| 2 dx, 



t € 



(3.12) 



(ii) Let A S H 1 ^). Then the solution above satisfies A £ C(R; H 1 ^)) and 

[\VA(x, i)| 2 dx + aln(l + |A(x,i)| 2 )] dx 

= / [|VA (x)| 2 dx + aln(l + |^ (x)| 2 )] dx, t G M(3.13) 

Proof. The norm conservations (|3.12p and (|3 . 1 3[) result from multiplying (|3.11[) 
by A and dtA respectively and integrating the complex and real parts respec- 
tively. This formal proof is justified by the standard truncation process. 

Let S(t) be the group operator associated to the linear Schrodinger equation 
idtA + AA = 0. Then the Duhamel formula for p. lip reads 



A(x,t) =SQk)Mx) 



(3.14) 



Since x i— > xj (1+x) is Lipschitz, we easily infer that the R.H.S. of (|3 . 14[) defines 
a contraction on a suitable ball of C([0,T]; L 2 (R d j) for some T > 0. The local 
well-posedness in L 2 (M. d )) follows. Global well-posedness is derived from the 
conservation law (|3.12|) . 

The H 1 theory follows the same argument, noticing that 



|A| S 



1 + |AP 



-A 



A 2 



(1 + \A\ 2 ) 2 



VA- 



\A\< 



(1 + \A\ 2 ) 2 



VA 



< - 2 m. 



□ 

Remark. As a consequence of ((3~T2l . (|3TT5|) and ln(l + \A\ 2 ) < \A\ 2 , we obtain 
the uniform bound 



|Vyl(x,i)| 2 dx < / |A)(x)| 2 dx 



|VA (x)| 2 dx, 



t £ 



(3.15) 



Contrarily to the context of the usual nonlinear cubic Schrodinger equation, this 
bound does not depend on the sign of a and in particular saturation means that 
no blow-up is occurs. 



3.2 Solitary waves — one-dimensional results 

In the one-dimensional case, it is possible to compute first integral formulations 
of the solitary waves. 

Bright solitary waves are sought for in the form A(x, t) = e %UJt u{x) (see [TT|). 
where A is a solution to (|2. 10|) with A^ = 0. The function u is supposed to 
have a maximum at x = (w(0) = u rn > and u'(0) = 0), therefore 

[u'(x)} 2 = (« - a)[u 2 {x) - u 2 J + a[ln(l + u 2 ) - ln(l + u 2 m )}. 



s 



We furthermore want that for x — ► oo, u(x) — > and u'(x) — » 0. This yields a 
unique possible frequency for the solitary wave, namely 

, , ln(l + Q 
tjj = a 1 ^ 



and 



Since w m is supposed to be the maximum of u, this quantity is positive only 
if a = 1 (focusing case) and the bright soliton is solution to the first order 
equation: 



[u'{x)f = a ( -^M ln(l + t&) + ln(l + u 2 ) 



u'ix) = -sign(z) Jln(l + w 2 ) - ln(l + v?J with w = 1 - ln ^+"^ , 
V ui 



Dark solitary waves are sought for in the form A(x,t) — u(x) (see [12] ) 
where A is solution to (12.10|) . There is no time-dependence. We assume that 
linXr—^oo u'(x) = and consistently with =/= 0, 

lim u[x) = — lim u(x) — Uoo. 

x — > + oo x — > — oo 

Then 

l + u 2 



At the origin u(0) = and we want more generally that |m(x)| < |woo|- There- 
fore, dark solitary waves only exist if a = — 1 (defocusing case) . In this context 
u(x) is a monotonous function and is solution to the first order equation: 



u'(x) = sign(u oc )y'u 2 -u^ - + In ( j^~2 
For both bright and dark solitary waves, no explicit solution is known. 



3.3 Solitary waves — a priori estimates and non existence 

Consider now the solitary wave solutions of (|3. 1 1|) in any dimension d, that is 
solutions of the type A(x, t) = e tuJt U(x), where U S H 1 (M d ) (we thus are only 
concerned with "bright" solitary waves). A solitary wave is a solution of the 
elliptic equation 

1771277 

-AU + wU = a±±^, UeH\R d ). (3.16) 

A trivial solution is U = 0. We seek other nontrivial solutions. 



9 



Lemma 2 Any H (Mr) solitary wave satisfies 



\U\ 



2 



VU\ 2 + a; -a . ' ,, |17| 



dx = (3.17) 



l + \U\\ 
(energy identity) 

(d-2) / \X7U\ 2 dx+duj [ \U\ 2 dx-ad [ [\U\ 2 - ln(l + |C/| 2 )] dx = (3.18) 

JRd JR 11 JB d 

(Pohozaev identity). 

Proof. As for Theorem [TJ (|3.17p results from multiplying (|3 . 1 6[) by 17 and inte- 
grating. To get ()3.18|) . one multiplies (|3.16|) by Xkd Xk Uk, integrates the real part, 
and sums from 1 to d. This is justified by a standard truncation argument. □ 

Corollary 3 No nontrivial solitary wave (solution of (|3.16[) ) exists when 

(i) a = — 1 (defocusing case), for lu > 0. 

(ii) a = 1 (focusing case) and to > 1 . 

(iii) a = ±1 i/w < provided \U\ 2 /(1 + \U\ 2 ) = 0(l/|x| 1+£ ) ; e > as |x| -> 



Proof. Identity (|3.17|) implies that no solitary wave may exist when a = — 1 
and w > or a = 1 and w > 1. When d — 1,2, Eq. (|3.18[) implies that no 
solitary wave exist when lu < and a = 1. Recall d = 2 is the physical case. 
The remaining cases (lu < 0, a = —1 or a = 1, d > 3) follow from the classical 
result of Kato[8] on the absence of embedded eigenvalues. Indeed, we can write 
(j3~T6l) as 

\U\ 2 



AU+ (-lu + V(tc))U = 0, l/(x)=a 



l + |i7| 2 ' 



assuming furthermore that V(x) = 0(l/|x| 1+e ), e > 0, as |x| — * oo. A proof 
for d = 3, 4 or d > 5, lu < — (d — 2)/2 with no decaying assumption is given in 
Appendix. □ 

Corollary 4 Solitary waves may exist only when a — 1 and < lu < 1 . 

Corollary 0] is consistent with the one-dimensional "explicit" result. We first 
have a classical regularity and decay result. 

Proposition 5 Let a — I and < to < 1. TTien any ?7 € F x (M d ) solution of 
(|3.16| satisfies 

UeH°°(R d ), (3.19) 
e 5 ' x '[7 e L°°(IR' i ) /or any (5 < w/2. (3.20) 
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Proof. U 6 H oa (R d ) results trivially from a bootstrapping argument using 
|£/| 2 /(1 + |£/| 2 )<1. To prove (pT20]) . we first derive the estimate 



/ e" |x| [\VU\ 2 + \U\ 2 ] rfx < +00. 

Js. d 



(3.21) 



In fact, as in CazenavejUIS], we multiply (|3 . by e"' x 'i7 and integrate the real 
part (this formal argument is made rigorous by replacing e"' x ' by e"l x l^ 1+e l x ", 
e > 0, e — > 0) to get 

/ e w|x| [\VU\ 2 + lu\U\ 2 ] dx< / e^M-Ml—dx + u; / e"l x l|[/||VC/|dx. 

(3.22) 

By (fXT5|) there exists i? > such that \U\ 2 /{1 + \U\ 2 ) < lu/4 on R d \ B R . Thus 
we infer from (I3.22p that 

/ e-l x l[|VC/| 2 +c|C/| 2 ]dx < / e«wJ^«fc+2/ e-l x l|C/| 2 dx 

+%(f e^\U\ 2 dx+ f e"l x l|V[/| 2 dx 

which implies (|3.2ip . 

Now we write U as a convolution 



\U\ 2 U 



U(x)=H u * ,\\ , where ff„ = T~ x ( ; |H2 ). (3.23) 



As it is well known ( UJ), ff w (x) = o)( d ~ 2 )/ 2 Gi(w 1/2 x) where 

'c\z\( 2 - d V 2 K (d _ 2)/2 (\z\), d>3, 
K Q {\z\), d=2, 



where K v is the modified Bessel function of order v. Furthermore (see pQ), one 
has the asymptotic behavior: 

K u (z) ~ \T{v){\\z\y v , for v > 0, as -»■ 0, 
tf (*) - -M|z|), as|z|^0, (3.24) 

if„(z) - C|^|-V2 e -NI 5 as|z|^oo. 

We infer from (|3~23)> that 

13 



e 5|x| |C/(x)|< / e 5 l x - x 'li/ w (x-x')e 5 l x 'l-^l— (x')dx' 

7R<i 1 + |C/| 2 



(3.25) 



Since by (13T241 e^Wii^x) G i 2 (M d ) for < <5 < w 1 / 2 , and by (I3T2H e^l |L/| 3 /(1-1 
|J7| 2 ) G L 2 (R d ) forJ < tu/2 < oj 1 / 2 , we deduce from ([3225]) that e s ^U G L°°(K d ) 
for < 5 < u/2. □ 
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Remark. Actually, the saturated cubic NLS equation should involve a small 
parameter e > 0, namely, in the focusing case, we should consider instead of 
IpTTTjl 

\A £ \ 2 A £ 



id t A £ + AA £ = 



(3.26) 



l + e|A e | 2 " 

Theorem Q] is of course still valid for a fixed e > 0, but (]3 . 13[) and (|3 . 1 5|) should 
be replaced by 



|V^l(x,i)| 2 dx + -iln(l + |A(x, t)\ 



rfx 



|VA (x)| 2 dx+^ln(l + |A (x) 



dx, 



f \X7A{x,t)\ 2 dx< - f \A Q (x)\ 2 dx + f \X7A {x)\ 2 dx. 

JR d e JR d JR d 

For solitary waves A £ (x,t) = e luJt U(x), (|3 . 26[) reduces to the elliptic equation 

\U\ 2 U 



-AU + ojU 



l + e\U\ : 



Setting V = e 1 ^ 2 !!, one obtains 



1 \V\ 2 V 
-AV + ujV=- 11 



e 1+ \V 



2 ' 



The only possible range for the existence of nontrivial H 1 solitary waves is 
u> G]0, l/e[. Proposition [5] is still valid for li in this range. 

3.4 Solitary waves — existence results 

We now turn to the existence of non-trivial H 2 solutions of 

\U\ 2 U 



-AU + uU = 



l + |f/P 



when < u> < 1. We will look for real radial solutions f(x) = n(|x|) = u(r) 
and thus consider the ODE problem 



// d - 1 / , 

l U + LOU 



1 + u 2 ' 



u e -ff 2 (]0,oo[), u'(0) = 0. 

We recall a classical result of Berestycki et al. [2] 



(3.27) 



Theorem 6 ( |2j, p. 143) Let gf 6e a locally Lipschitz continuous function on 
M + = [0, +oo[ such that g(0) = 0, satisfying the following hypotheses. 



(HI) a = inf{C > 0, g(Q > 0} exists and a > 0. 
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(H2) Let G(t) = f*g{s)ds. There exists ( > such that G(() > 0. 

Let Co — inf {C > 0, ^(C) > 0}. Ln view of (HI) and (H2), Co exists and 
Co > a- 

(H3) lim s ^ a g(s)/{s-a) > 

(H4) g(s) >Oforse]a,Co]. 

Let (3 = inf{C > Co, G(C) > 0}. Jn mew o/ (H4), Co < < +oo. 

(H5) If [3 = +oo, tften lim 5 ^ +00 9 (s)/s* = /or some Z < (d + 2)/(d - 2) f«/ 
d = 2, we may choose for I any finite real number). 

Let us consider the Cauchy problem 



Then there exists C G]Co,/3[ such that (]3 . 28[) has a unique solution satisfying 
u(r) > for r e]0, +oo[, u'(r) < for r g]0, +oo[ and lim,.^ +00 u{r) = 0. If 
in addition limsup^g g(s)/s < 0, then there exists C > and S > swc/i i/iai 
< u{r) < Ce~ Sr , for0<r < +oo. 

Theorem 7 If a = 1 and < u> < 1, £/iere exists a nontrivial positive solution 



Proof. The case d = 1 has been addressed in Sec. 13.21 Consider now d > 2. 
We apply Theorem [5] with 



which graph is displayed in Fig. [T] 

Note that a = \Juo/{l — w) which yields (HI). Setting u = y/uj/Jl — oj) +e, 
one easily checks that 




(3.28) 



of (£23). 



s(") 



2w(l - w) + O(e), 



it — a 




a 



Figure 1: Graph of function 5 
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and (H3) is satisfied. One computes G(u) = (1 — u)u 2 — i ln(l + u 2 ), which 
obviously satisfies (H2) and (H4) with [3 = +00. Last (H5) holds true (for 
I > 1). □ 

Remark, u satisfies the decay rate of Proposition \5\ 

4 The Zozulya— Anderson System 
4.1 Estimate on the potential 

We now restrict to the space-dimension d = 2 which is the context of the 
derivation. To mimic the proof for the Cauchy problem in the one-dimensional 
case, we would like to express ip in terms of A for say A 6 L 2 (K 2 ). With such a 
data A, we indeed have a unique ip in some convenient space but no Lipschitz 
regularity for the mapping A 1— ► cp, which is required to perform some fixed 
point procedure. To ensure this we will have to assume A G H 2 (R 2 ). 

To derive the first estimates, we consider time as a parameter and do not 
express it. We therefore introduce the weighted homogeneous Sobolev space 

H = {ip e S'(R d ), (1 + \A\ 2 y/ 2 Vip g L 2 (R d )}/R 

together with its natural Hilbertian structure. 

Lemma 8 (i) Let A G L 2 (M. 2 ). There exists a unique tp G H solution of 

div((l + \A\ 2 )Vp) = d x (\A\ 2 ) in V'(R 2 ) (4.29) 

such that 

[ {l + \\A\ 2 )\Vp\ 2 d*< 1 - f |A| 2 dx. (4.30) 

(ii) If furthermore A G H 2 (W 2 ), then Vp G H 2 (M. 2 ) and there exists a polyno- 
mial P vanishing at such that 

||V^|| ff 2 (R2) < P(\\A\\ H 2 m ). (4.31) 

Proof, (i) We define a smoothing sequence (9 e ) e> o with J R2 9 £ dx — 1 and 
A e = A * 6 e is such that A e — > A G L 2 (R 2 ). In particular 

||^4e||L 2 (R 2 ) < II^IIl 2 (r 2 )- (4.32) 
By Riesz theorem there exists a unique solution to 

div((l + |A £ | 2 )V^)=^(|A e | 2 ), (4.33) 

i.e. 

- Ap e - div(|A £ | 2 V^) = d x (\A £ \ 2 ) (4.34) 

after noticing that the R.H.S. of Eq. (|4.33[) defines a linear continuous form on 
H given by 

(d x {\A £ \ 2 )^) = [ \A E \ 2 d x ^dx. 

Jm 2 
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Now we get from (|4.33[) 

/ (l + |A e | 2 )|V^| 2 dx = - / \A E \ 2 d x ^dx 

JR 2 JR 2 

< \ f \A £ \ 2 dK+ \ [ \A £ \ 2 (d x <p E ) 2 dx, 

1 JR 2 1 JR 2 

which yields (together with (|4.32[0 

/ (l + i|A e | 2 )|V^| 2 dx<i / \A £ \ 2 dx<U \A\ 2 dx. (4.35) 

JR 2 * z JR 2 z JR 2 

Up to the extraction of a sub-sequence, we have V<^ e — > and d x (|A e | 2 ) — * 
d x (\A\ 2 ) in V'(R 2 ). From Eq. (|4T3"5]) . A e V<^ £ B weakly in L 2 and for all 

/ A e W(p e ■ Wipdx -> / A\Jip ■ Vipdx, 

JR 2 JR 2 

therefore B = AVip. Since ||A e | 2 V(^ £ | = | vl e j j^4 e Vy? e ] , \A £ \ 2 Vcp e -> |A| 2 V<^ in 
I?'. We can pass to the limit in Eq. ()4.34[) and obtain 

-Atp - div(|A| 2 V</?) = ^(lAI 2 ) in P'(K 2 ), 

i.e. div((l + \A\ 2 )Vtp) = d x (\A\ 2 ) and deduce estimate (!4~30l) from (|4~35l This 
yields the existence of ip £ H. The uniqueness is straightforward: two solutions 
tpi and v?2 would satisfy 

/ (1 + i|A| 2 )|V(^i - <^ 2 )| 2 dx = 0, i.e. V((^i - <&) = a.e. 

and hence be equal in i?. 

(ii) We first notice that |j4| 2 A</? is meaningful in H^ 1 (R 2 ). Actually, for any 
ip € H 1 ^ 2 ), one defines 

(|A| 2 A(p,l/')H-i(R 2 ),-ff 1 (K 2 ) = (Ap^Al 2 ^) H-i(R 2 ),m(R 2 ), 

which makes sense since \A\ 2 ip £ iJ^R 2 ) for A € H 2 (R 2 ), V € i/^M 2 ). Thus 
we can write (|4.29[) as 

(1 + \A\ 2 )A^ = -V\A\ 2 ■ + d x (\A\ 2 ), 

and 

We claim that F € L r (R 2 ), for any r € (1, 2), with 

||^1|l>'(K 2 ) < C||A|| i 2( R 2)|| J 4|j fl -2( R 2). 

First, |V|^| 2 • V<^| < 2|VA||AV^| and by Holder 

IIVjAI 2 • \/ip\\ L r m < 2||VA|| i p (R3) ||AV^|| L 2 (K 2 ) 
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for any 1 < r < 2 and p = 2r/(2-r) e (2,oo). Since || AV( / 9|| L 2 (R 2) < ||A|| L 2 (R 2) 
and tf^IR 2 ) C L«(M 2 ) for all g > 2, we obtain that 



V|4| ; 



1 + \A\> 



< c||^IIh 2 (r 2 )II^IIl 2 (r 2 ), 



1< r < 2. 



Similarly 



1 + |4| 2 



L r '(R 2 ) 



< 2||A0 x 4|| L r (R 2) < 2||i4|| L 2 (R 2 ) ||0 I A|| i p (R 2 ) 



< C||A|| H 2 (K 2 ) ||A|| L 2 (R 2 ) , Kr<2. 
By elliptic regularity, we infer thus that for any r, 1 < r < 2, 

||V<^|| W l,r( K 2) < C||A|j fl -2( R 2)|jA|j i 2( R 2). 

By Sobolev embedding, 

||V^|| L9 (R2 ) < C||V^|| W l,r(R2) < C||4|| H 2 (R 2)||A|| L 2 (R 2). 

for i = i - i i.e. g = 2r/(2 - r) for all r, 1 < r < 2. Thus for any p > 2 



g r 2 ' 

V|4| 2 



1+L4l; 



< 



LP(R 2 ) 



l|V^|| L 2 p(R 2 ) |)V|^| 2 || L 2p (R 2 ) 



< C||A|| ff 2( R 2) 11^411^2^2) 11-411^2^2) 



C\\A\\ LHR2) \\A\\ 

H 2 (R 2 



(we have used the fact that H (K ) is an algebra and the embedding H 



C 



2 ) for allg>2). 
Similarly, for any p > 2 

5,(|4| 2 ) 



1 + L4f 



< 



Lp(R 2 ) 



2||4|| L 2 P(R2) ||a x A|| L 2p 



< 



cii^4iifj' i (R 2 )imiiJ 2 (M 2 )- 



Finally for any p > 2 

||F|| iP(K 2) < ^11411^(^,(1 + 11^^2^2)1141^2^2)). 

and by elliptic regularity 

||V^|| W l,P (R 2) < C||A|| 2 f2(K 2)(l + m|| L 2 (R 2 ) ||A|| Jf 2 (M 2 )), Vp > 2. 



We now check that Vtp ■ V|4| 2 /(l + \A\ 2 ) e H 1 ^ 2 ). This easily reduces to 
showing that V(V<£ ■ V|4| 2 ) e L 2 (R 2 ). For e {1,2}, 9 Xi 9 x ^ e L P (M 2 ) 

since 

?2 ), P 



d Xi d x (f = T^rfA^ and A</j e 



Id 2 
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Thus d Xi d Xj (pV\A\ 2 G L 2 (R) (V|A| 2 G H^R 2 ) c L«(R 2 ), Vg > 2). 

On the other hand, taking p > 2 we see that V<p G L°°(R 2 ) and thus 
V^AJAI 2 G L 2 (R 2 ). 

It is also easy to check that d x {\A\ 2 )/(l + \A\ 2 ) G PT^R 2 ). 
Finally, Atp = F G i/ 1 (IR 2 ), proving that Vy> G P 2 (R 2 ) with an estimate of the 
form 

||V<£>||iJ2(R2) < P(m|| ff 2( R 2)), 

where P is a polynomial vanishing at 0, which proves (|4.31[) . □ 

Remark. All above estimates are therefore uniform in time, and if A G C([0, T]; H 2 (M. 2 )) 
for some T > 0, one has 

l|V^||c([0,T];_f/ 2 (R 2 )) < P(\\A\\c(lO,T];H 2 (R 2 )))- 

4.2 Solitary waves — non existence results 

We now look for solitary wave solutions of (|2.9p , that is solutions of the form 
(e iut U(x), 4>{x)) with x G R d , lu G R, U G P 1 ^), and $ € H. Thus ([/» 
should satisfy the system 

f -At/ + wP = aUd x d>, 

I (4.36) 
\ div((l + |[/| 2 )V^)=^(|[/| 2 ). 

The existence of nontrivial solutions of (|4.36p is an open problem. Note that 
(|4.36p does not seem to be the Euler-Lagrange equation associated to a varia- 
tional problem. We have however: 

Proposition 9 (i) Let a = — 1 (defocusing case). Then no nontrivial solution 
of (|4.36[) exists for lu > 0. 

(ii) Let a = 1 (focusing case). No nontrivial solution of (|4.36[) exists for lu > 1. 
(hi) Let a = ±1. No nontrivial solution of (|4.36[) exists if lu < provided 
d x( f> = 0(l/|x| 1+£ ) ; e > as |x| -> +oo. 

Proof. From (|4.36p we have 

/ \VU\ 2 dx + Lu { \U\ 2 dx = a [ \U\ 2 d x <pdx, 

JS. d J& d J& d 

(i + |c/| 2 )|V0| 2 dx = - f \u\ 2 d x( f>dx, 

jR d JR d 

and 

f \VU\ 2 dx + iu [ \U\ 2 dx-a[ (1 + |P| 2 )|V0| 2 dx = 0, (4.37) 

J«. d JM d JS. d 

which proves (i). Now independent of the sign of a, and from (|4.30p and (|4.3T[) . 

/ \VU\ 2 dx + Lu [ \U\ 2 dx< [ \U\ 2 dx. 
Jm d Js. d JR d 

Thus 

{ |VP| 2 dx + (lu- 1) ( |P| 2 dx<0, 

JK d JM d 

which proves (ii). Part (iii) results from [5]. □ 
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4.3 The Cauchy problem 

We consider the system 

id t A + AA = -aAd x <p, 

div((l + |A| 2 )V^) - ^(|A| 2 ), (4.38) 

A(;0)=A . 

Theorem 10 Let A a E H 2 (R 2 ). 

Then there exists T > and a unique solution (A, Vy) of (|4.38|) such that 
A e C([0,T ];iJ 2 (R 2 )) andVip eC([0,T ];H 2 (R 2 )). Moreover 

\\A(-,t)\\ LHM2) = ||Ao|| L 2 (R2) , 0<t<T a 

and 

f (1 + i|A| 2 )|V^| 2 dx < \ [ \A \ 2 dx, 0<t<T o . 
Jm. 2 z 1 Jm 2 

Proof. Uniqueness. Let (A, V<^) € L°°(0, T; iJ 2 (lR 2 )) and (B, Vtp) G L°°(0, T; ff 2 (R 2 )) 
two solutions of (|438"|) with 0) = B(-, 0). Then from fl225]) 2 one gets 

A(^ - v) + div(|A| 2 v^ - |s| 2 w.) = a,(|A| 2 ) - a x (|i?| 2 ), 

yielding 

|V(p- V)| 2 rfx+ / |A| 2 |V(^- V)| 2 dx 



(IA^-ISH^^-^x- / (\A\ 2 -\B\ 2 )^-V^~^)d> 



(4.39) 

Observing that |A| 2 - \B\ 2 = A(A - B) + B(A - B), the R.H.S. of (|439l is 
majorized by 

- ( \d x ^-^)\ 2 dx+{\\A\\ L ~ {m + \\B\\ L ~ m ) [ \A-B\ 2 d x 
4 J& 2 Jr 2 

+ - j |V(^-^)| 2 dx 

+ (||A|| i « B» ) + ||B|| £ op (Ra) )||V^||Lo. (R3 ) / \A-B\ 2 dx 

Jr 2 

and by Sobolev embedding 

l|v(v? - V0IU 2 (m 2 ) < cdl^lU^E 2 ), ||^IIh2 (R 2), ii v^ii^^))!!^ - bw^^. 

(4.40) 

On the other hand, we obtain readily from (|4.38jl i that 



ld_ 

2 dt 



f \A-B\ 2 dx< f \A-B\ 2 \d x <p\dx+ f \B\\B x (<p-il;)\\A-B\dx 
Jm 2 Jr 2 Jr 2 
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which together with (|4.40p and the Cauchy-Schwarz lemma yields 

< C(l!^ll_f/ 2 (R 2 )j ll^llff 2 (R 2 ), ||V^||ff2( R 2), |jV-0||_H-2( R 2))||A - B|| £ 2( R 2) 

and A = B by Gronwall lemma. 

H 2 a priori estimate. We derive a (formal) H 2 a priori estimate on the solution 
of (|4.38| . Since _ff 2 (M 2 ) is an algebra, we deduce from Lemma [8] that 

\\^9a;ip\\c([0,T];H 2 (M 2 )) < ll^llc([0,T];i/ 2 (E 2 ))-P(ll^llc([0,T];i/ 2 (R 2 ))), ( 4 -41) 

where P was introduced in (|4.31|) . From the energy estimate 

~\\M;t)\\ 2 H *(R 2 ) < C\\AdM-MH><F>)\\A(-,t)\\H*(*»), 
we infer with (|4.41|) the local H 2 bound 

\\A(;t)\\ H 2 m <C(\\Ao\\ H 2 m ) forO<t<T , (4.42) 
T < T sufficiently small. 

Approximation of (|4.38[) . The strategy is now to implement a compactness 
method using the (justified) a priori estimate (|4.42j) . For e > 0, we consider 
the system 

id t A £ + AA £ = -aA £ d x <p £ , (4.43) 
div((l + £ A 2 + |A £ | 2 )V^) = d x (\A £ \ 2 ), (4.44) 
A £ (-,0) = A . (4.45) 
Solving V<p e in terms of A e , we obtain from (|4.44[) that Vip £ satisfies 

ef |AV^ E | 2 dx+ / (l + \\A e \ 2 )\Vtp £ \ 2 d*<\ ( \A £ \ 2 dx. (4.46) 
Jr 2 Jr 2 2 2 J R 2 

Well-posealness of approximate system. We now check that the Cauchy problem 
(j4^3l) - (|4^431) is globally well-posed in H 2 (R 2 ). Let first A £ , B £ e H 2 {R 2 ) and 
(p £ , i/> £ the corresponding solutions of (|4.44|) . Proceeding as in the uniqueness 
proof above, one gets 

ef |VA(^ E -ip e )\ 2 dx + f |V(^ £ -V e )| 2 rfx 
Jr 2 Jr 2 

< C(II^ £ ||h 2 (R 2 )i \\B e \\H 2 (R 2 ), ||V^ £ ||h 2 (R 2 ), ||V^ e ||ff2( R 2))||yl e - i? e |j L 2( R 2). 

(4.47) 

Denoting d x ip £ by F £ (A £ ) we write (|4.43|) on the Duhamel form with S(t) — 
exp(iiA), 

A £ {t) = S(t)A -a [ S(t- s)A £ F £ (A £ )ds. (4.48) 
Jo 
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Using (|4.47j) and the unitarity of S(t) in iP(R 2 ), we deduce that the R.H.S. of 
p~4"8)) defines a contraction in C([0, T e ]; H 2 {R 2 )) for some T £ > 0. 

This implies the local well-posedness of (|4.43l) (|4.45[) in H 2 (R 2 ). Using the 
H 2 bound P~4"T)) on Vtp £ , we infer from (|4~4^|) an a priori bound in C([0, T]; H 2 (R 2 )) 
for A e and for all T > 0. This proves that the Cauchy problem (|4.43|) — (|4.45[) is 
globally well-posed, for any fixed e > 0. 

Limit s — ► 0. Now we have the bounds (I4.46P and 

||VvJ e (-,t)||z"(B») + Vi||AV^(.,t)||L2 ( B»)<C 1 0<i<T, (4.49) 

where C and T do not depend on e. Moreover, from (|4.43p and (14.49[) we have 
a bound on dtA e which is independent of e: 

||9 t A e (-,i)|| i2(R2) <C, 0<t<T. 

It is now standard to pass to the limit as e — > (see Q33). By the Aubin-Lions 
compactness lemma, we obtain a subsequence (A £ ,\7ip e ) such that A £ — » A 
in L°°(0,T;iJ 2 (R 2 )) weak-star and L 2 (0, T; fl^QR 2 )) strongly, V<^ £ -> Vvs in 
J L oo (0,T;iJ 2 (R 2 )) weak-star and i 2 ([0,T] x K 2 ) weakly. The limit (A,V<p) 
belongs to (L°°(0, T; _ff 2 (R 2 ))) 2 and satisfies (j438| . In fact (f4~38]) 1 is satisfied 
in i 2 (R 2 ) and (|4~38]) 9 is satisfied in iJ^R 2 ). 

The fact that (A,V<p) e (C(0,T; H 2 (R 2 ))) 2 results from the Bona-Smith 
approximation (see [3]). □ 

Remark. We do not know whether the local solution obtained in Theorem ITOl 
is global or not. 



5 Conclusion 

We have given a full description of how to derive from the Kukhtarev equations 
an asymptotic model for the propagation of light in a photorefractive medium. 
This derivation is only heuristic insofar as asymptotics are not justified, which 
would be out of reach now. Some properties of photorefractive media such as 
memory have also been neglected. 

The ID asymptotic model is a saturated nonlinear Schrodinger equation the 
Cauchy problem of which is studied (in any space dimension) in L 2 and H . 
We also prove the existence of solitary waves in one and higher dimensions. An 
interesting and open issue would be to study the transverse stability of the ID 
solitary waves in the framework of the asymptotic model. 

For the 2D asymptotic model (the Zozulya- Anderson model) we also have 
studied the Cauchy problem and the non-existence of solitary waves. The ques- 
tion of imposing other boundary conditions, not vanishing in one space direction, 
can also be addressed to treat a wider range of experimental applications. 
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A Non-Existence of Solitary Waves in Non-Physical 

Cases 

The goal is here to complete the results of Corollary [3] for u> < with no decaying 
assumption. We have already seen that Eq. (|3. 18|) implies that no solitary wave 
may exist for d = 1, 2 and a = 1 (focusing case). 

To go further, let us use both Eqs. (|3.17|) and (I3.18|) to obtain 

/ f2o,+ (d - 2 |y )\U\ 2 dx-ad[ [|lf -ln(l + |Lf)]dx = 0. 
We set 

f(x) = hu + ^y^x ) x ~ ad ^ x ~ ln(1 + 

and we know that f Rd F(|£/| 2 )dx = 0. Now F(0) = and 

2X 2 (u -a)+ X(Au - (4 - d)a) + 2uj 

F {x) (TTxT < ' 

if uj < 0, a = 1 and d — 3,4. Therefore F(|C/| 2 ) = a.e. By a bootstrapping 
argument, we notice that any H 1 solution to Eq. (|3.16[) is indeed in H k for all 
k and therefore continuous. Hence F(\U\ 2 ) = on R d . Since F'(X) < the 
only possible value for U is U — on M. d . 

We can refine this result, finding other parameter ranges for which 2X 2 (uj — 
a)+X(4uj-(4-d)a) + 2uj < 0. If d > 5 and a = 1, this holds for oj < -(d-4)/4. 
Moreover, 



2X 2 (w - a) + Jf (4a> - (4 - d)a) + 2w 
= 2 



(oj - a)(X - l) 2 + X(4lu - 4a + ^a) + w + a 



No solitary wave can exist for a = — 1 and u; < — 1. Hence we complete Corollary 
[3] with 

Corollary 11 No non-trivial solitary wave (solution of (|3.16|) ) of the saturated 
NLS equation exists when 

(i) a = — 1 (defocusing case), for u) < — 1. 

(ii) a = 1 (focusing case), for u> < 0, if d — 3, A and lu < — (d— 4)/4 if d> 5. 
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